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(~) Abstract. We give an update on the lifetime difference of B s mesons which accounts 

| for recent lattice bag parameter results, and obtain (Ar/r)s s = (9.3^$)%. We then 

■ review the current theoretical uncertainties and conclude with a pessimistic perspective 

Oh, on further improvements. 

1. Status of theory 

X. 

The width difference in the B s system is expected to be the largest rate difference 
in the B hadron sector [l|]. Some experimental bounds on this quantity from LEP 
@ and CDF || exist, and in the near future ATb s will be measured quite precisely 
@. Several factors contribute to the interest in Ar Bs : a large value of the width 
difference opens up the possibility for novel studies of CP violation without the need 
for tagging |4| . Moreover, an experimental value of Ar# s would give information about 
the mass difference in the B s system || (although at this moment it appears that the 
mass difference will be measured sooner than the width difference). Another interesting 
point is that new physics can only lead to a decrease of the width difference compared 
to the standard model value |J . An experimental number which is considerably smaller 
than the theoretical lower bound, would thus be a hint for new physics that affects 
B s -B s mixing. Besides the need for a reliable theoretical prediction of /ST S in order to 
fulfill the above physics program it is of conceptual interest to compare experiment and 
theory in order to test local quark-hadron duality, which is the underlying assumption 
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in calculating heavy quark decay rates. One can show that duality holds exactly in the 
limit Aqcd *C m b — 2m c <C m h and N c — > oo J7|. So far no deviation from duality 
has been conclusively demonstrated experimentally and theoretical models of duality 
violation in B decays tend to predict rather small effects [|J. 

During the last year there has been remarkable interest in the lifetime difference 
of B s mesons. Besides several experimental studies, many lattice calculations of the 
relevant nonperturbative constants were done |§, [|K|, ||11|| . This improvement in theory 



input motivates the present update of the result presented in |12|]. We also clarify the 
origin of seemingly disagreeing recent evaluations of Ar^ s . 

The theoretical status of the lifetime difference computation is as follows: the heavy 
quark expansion (HQE) allows us to expand the decay rate of heavy mesons in inverse 
powers of the heavy quark mass. The leading term is described by the decay of a free 
quark (the so-called parton model) and therefore equal for all B hadrons. The next term 
is suppressed by two powers of the heavy quark mass; it is related to the kinetic and 
the chromomagnetic operator. The so-called weak annihilation and Pauli interference 
diagrams contribute first at third order. Operators involving the spectator quark begin 
to appear at this order and these are mainly responsible for the lifetime differences of 
B hadrons. (A smaller contribution comes from the matrix elements of the kinetic and 
chromomagnetic operator; for the lifetime difference of B s mesons these terms vanish 
identically.) We can write for the decay rate difference of two mesons 



A 3 

Ar = — 



r ^^ 1, + -) + ^ 0, + -+) + 



4/T ' / 

The IYs are products of perturbatively calculable Wilson coefficients and matrix 
elements, which have to be determined by some non-perturbative methods like lattice- 
QCD or sum rules. For a comparison of experiment and theory we need besides the 
experimental value several theoretical ingredients. First, the perturbative prediction in 
leading order (r^), second, corrections due to the strong interaction (rjj : NLO-QCD), 
third, subleading l/m b -corrections (I4 ) and last, but not least, the determination of 
the appearing non-perturbative parameters. For (AT/T)b s we have all these pieces: the 
leading term was calculated in |13[ in the factorization approximation; 1 / m^-corrections 
were determined in [14]; the ct s -corrections are given in [12]; and the lattice values for the 
decay constant fs s and the bag parameters B can be taken from computations directed 
at the mass difference. Bs, which is specific to the width difference was calculated in 



10 1 , [11]. (The matrix elements of the dimension- 7-operators, which emerge in r 4 , 



have been estimated up to now only in vacuum insertion approximation |14]|-) As the 
calculation of was the first to consider QCD corrections to spectator effects, there 
is also a conceptual point of interest. Soft gluon emission from the spectator s quark 
leads to power-like infrared divergences, which would spoil the HQE. In [I2| the infrared 
safety of Ar s was explicitly shown, in agreement with the theoretical argument of | |15| |. 

Compared to Ar s our knowledge about the lifetimes ratios T(B + )/r(Bd) and 
T(Afc)/r(I?d) is much poorer. Here only I3 is known |L£|, while the calculation of 



and is still missing. Moreover we only have preliminary lattice studies for the 
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bag parameters in T 3 Results from sum rules for the bag parameters are discussed 
in 1113. For these lifetime ratios more work needs to be done. 



2. Numerical update and uncertainties 



The following update is based on the NLO expressions given in |TI|. The width difference 
is normalized to the semi-leptonic B s branching fraction; numerically we find 



230 MeV 



0. 007 B(m b ) + 0.132 



M 2 B B s (m b ) 



-0.078 



(m b + m s ) 2 

We factored out the decay constant and use B and Bs to parametrize the matrix 
elements 



(B s \(b lSl ) V - A (b ]Sj ) v „ A \B s } = ~f B M 2 B B, 



(5 s |(6 i Sj) s _p(6 i s i ) s _ P |S s ) 



Ml 



:B 



s- 



(2) 
(3) 



3s (m b + m s) - 

- 4.8 GeV. The third term in 



Here rh q denote MS quark masses at the scale m b 
square brackets in Eq. (|1|) is an estimate of the l/m b correction in the factorization 
approximation [14]]. (This correction is slightly larger than in [12|, because we now use 
fh s = 0.1 GeV. We will also use fh b = 4.2 GeV.) We note that 

i) the term involving the parameter B(m b ), which also appears in the mass difference, 
is negligible; 

ii) the NLO correction to the coefficient is large and reduces the width difference. 
In LO we obtain the coefficients 0.011 and 0.203 instead of 0.007 and 0.132, 
respectively; 

iii) the l/m b correction is also large and negative, and its importance is amplified by 
the negative NLO correction. 

As a consequence recent estimates of the width difference tend to be significantly 
smaller than the leading order estimate of [14]] and those based on the factorization 
approximation (B = Bs = 1). 

There is another way of representing the result ([!]) [|TU[ based on the observation 
that B{m b ) appears in the mass difference an the fact the mass difference for B^ mesons 
is accurately determined experimentally. Then 



\t b Am B , 



M, 



B A 



Vt 



Vt 



td 



0.030 - 0.937R s (m b ) - 



0.35 
~B~ 



•(4) 



In this representation the value of the first bracket can be taken from experiments, and 



K 
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V cb V cs 



v ts v tb 



1 



T] B (m b )So(x t ) ^ 
is a known factor, if we assume that the CKM matrix is unitary. The advantage of this 
representation is that hadronic uncertainties enter only in ratios, i.e. in 



fB s B B s ( m b) 
fL B B d ( m b) 



R s {m b ) 
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B s {m b ) 



8 (m b + m s ) 2 B{m b ) 
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B{m b ) 


Bs(m b ) 


fB s /MeV 


(Ar/r) Bs 


BBGLN98 [12 


0.90 


1.07 


210 


6.0% 


Hashimoto et a/. || 


0.85 ±0.11 


1.24 ±0.16 


245 


10.7% 


Becirevic ei a/. ]TU| 


0.91(3)18:88 


1.32(3)18:8! 


uses Eq. (g) 


4.7% 


Gimenez/Reyes |IH 


0.83 ±0.08 


1.25 ±0.14 


uses Eq. (g) 


5.1% 


this summary 


0.9 ±0.1 


1.25 ±0.1 


230 ± 10% 


9.3% 



Table 1. Summary of recent evaluations of (Ar/r)s,. For theoretical uncertainties 
on (Ar/r) Bs consult text. B s (m b ) = M| B s (m b )/(rh b + fh s ) 2 . 



and these are believed to be better known than fs s and Bs(m b ). This advantage 
is more than compensated by the need to know |Vf S /14d| in Eq. (||), which makes 
the prediction for ATb s sensitive to the global fits to the unitarity triangle and the 
theoretical assumptions that go into it. In particular, the prediction for AT Bg now 
depends on the assumption that the standard model describes flavour mixing correctly, 
even though the decay of B s meson is unlikely to be affected by new physics in mixing. 
For this reason we would rather discourage the use of this method to obtain Ar_B a . 

In Table [I] we present a summary of recent estimates for ATb s , based on either 
Eq. fl]]) (default) or Eq. (g). (We use here only estimates based on lattice calculations of 
B(m b ) and Bs(m b ).) It is evident that there is remarkable agreement on the hadronic 
parameters B(m b ) and Bs(m b ), and that the spread of results is mainly caused by 
different values of fs s or the use of Eq. (f|). In the following we first assess the theoretical 
uncertainties of the NLO calculation and then present our best evaluation, given as "this 
summary" in Table [l]. The current theoretical uncertainties are as follows: 

i) Residual scale dependence. This dependence arises from residual effects in the 
matching of the NLO Wilson coefficients for Ar^ to the parameters B and Bs 
computed on the lattice, and from the matching of the NLO Wilson coefficients for 
Are s to the Wilson coefficients of the AB = 1 weak effective Hamiltonian. Here 
we estimate only the second source of scale dependence, using the coefficients for 
scales m b /2 and 2m b {m b = 4.8 GeV) given in [0. We then find 

V T A calc [-2.7%) 1^230 MeVj (m b + m s Y ' {) 

ii) Normalization. There is an overall normalization error in Eq. ([I]), which results 
from the uncertainty in the measurement of the semi-leptonic branching fraction 
and from the value of the b quark mass. We estimate this error to be 10%. 

iii) l/m b - correction. The second major source of theoretical uncertainty originates 
from the term "—0.078" in Eq. (g), which has been obtained using the factorization 
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this uncertainty, we introduce four parameters that measure the deviation from the 
factorization approximation, corresponding to the four independent operators at 
order l/m b . We find that all l/m b effects add essentially constructively, so that 
there is no reason to suspect an amplification of corrections to the factorization 
approximation as a result of cancellations in the factorized expression. We then 
find 

where Ar parametrizes a typical deviation from the factorization approximation, 
corresponding to quantities like B — 1, and Bs — 1 in leading power. There exists 
no reliable information on violations of the factorization assumption for the matrix 
elements in question, but experience with leading order matrix elements suggests 
that Ar « 0.3 is a reasonable upper limit. 

iv) B{m b ) and Bs{m b ). Substantial work has been invested in controlling the nonper- 
turbative parameters at leading power. Table |I] shows that the results are rather 
consistent with each other (the number 1.07 in the first line is obsolete) and suggests 
that 

/ x - , s Ml B s (m b ) , , 

B(m h ) = 0.9 ±0.1, B s (m b )= -. Bs M 1.25 ± 0.10. 9 

(m b + m s y 

We note that at the present stage the uncertainty in B and Bs has become a 
minor factor, the dominant ones coming from residual scale dependence and l/m b 
corrections. 



We can now combine these uncertainties with Eq. (P to obtain our final result 

(t-)b. = (^^)\ 1±0 -0[( 13 - 2 --) % ^^) + °- 7%i? ^) 

- (7.8 ±1.8%)" 

= (9-3^)%. (10) 

The final number is obtained by adding all errors in squares and assigning a 10% 
uncertainty to fs a = 230 MeV. It is worth noting that this number is considerably 



larger than the result obtained in [|10| , |TT| for two reasons: first, the overall normalization 
obtained from Eq. ([|) is about 25% smaller compared to that of Eq. (|l|) with fs a = 
230 MeV, since the global CKM fit prefers a smaller value of the B meson decay constant; 



second, Ref. [|10| uses m b = 4.6 GeV, which increases the l/m b correction by almost a 
factor 1.5. 



3. Prospects for improvement 



Despite (or because of?) extensive work on radiative and l/m b corrections the theoreti- 
cal prediction of the width difference of the B s mass eigenstates remains rather uncertain. 
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This is due to an unfortunate conspiracy of negative corrections at next-to-leading order 
in a s and in the heavy quark expansion. The large error exhibited by Eq. fllUP makes 
it improbable that new physics would be first observed in the width difference, since 
a new physics contribution to the B s -B s mixing phase large enough to cause an effect 
larger than the theoretical uncertainty would rather been seen elsewhere, for instance 
as a time-dependent asymmetry in B s — ► J/i/jcj). 

Is further improvement possible? It is conceivable that future lattice calculations 
will reduce further the uncertainty in B and Bs, but the impact of the present 
uncertainty in these parameters on the total error is no longer dominant. Some 
improvement could be possible concerning the overall normalisation, but this effect 
can also not be substantial. The obvious targets for improvement are therefore 

ii) scale dependence: its reduction would require the calculation or estimate of a? s 
corrections. This appears to be a hard endeavour, and it would also require the 
computation of three-loop anomalous dimensions of the weak effective Hamiltonian. 

ii) l/m& corrections: obviously, the relevant matrix elements should be computed 
on the lattice to rid the calculation from the factorization assumption. Such a 
calculation does not appear likely soon, given the familiar difficulties with higher 
dimension operators. Since the l/m& corrections are large, one may think of 
computing the a s corrections to them. While this is feasible, the effort is probably 
too large as long as the matrix elements of the operators are not known accurately 
enough. 

It therefore appears that (Ar/T) b s will remain considerably uncertain for the foreseeable 
future. 
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